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ABSTRACT 


Experimental  values  for  recombination  rates  between  electrons  and  cesium  ions 
are  presented  in  this  report.  Measurements  cover  temperatures  from  4300°K 
to  1200°K  at  electron  densities  from  1.  5  (10^)  to  2  (10iZ)  cc~*.  The  cesiurr. 
plasma  is  generated  by  a  high-velocity  steel  spherical  projectile  fired  into  a 
special  chemical  kinetics  ballistic  range.  Microwave  F-band  and  Ka-band 
focused  probe  measurements  of  electron  density  were  made  as  a  function  of 
time  at  two  wave  locations.  Measured  and  calculated  electron  density  decays 
are  shown  to  agree  satisfactorily.  The  experimental  recombination  rates  are 
compared  to  tl  eoretical  values  from  Bates,  Kingston,and  McWhirter;  agree¬ 
ment  is  reasonable. 


The  experimental  apparatus  and  techniques  are  described  in  detail  in  the  report. 
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I.  INTRODUCTION 


Interest  in  alkali -metal  ion  recombination  with  electrons  at  Avco  RAD  was 
stimulated  by  discrepancies  among  available  experimental  rate  constants. 
Previous  measurements  by  Mohler,  ^  Dandurand  and  Holt^,  Wada  and  Knecntli^ 
D'Angelo  and  Rynn“l,  and  King^,  are  summarized  in  table  I 

Measurements  from  references  1,  j ,  4,  5  were  all  carried  out  at  approximately 
2000° K  Although  the  electron  densities  for  this  group  of  measurements 
varied  by  a  factor  of  4000,  the  recombination  constants  were  of  the  same  order 
of  magnitude.  The  room  temperature  measurement  in  reference  2  is  unique; 
it  appears  to  belong  to  an  entirely  different  mechanism  than  the  higher  temper 
ature  values. 

The  problem  of  ion-electron  recombination  has  be  on  treated  theoretically  by 
Bates,  Kingston,  McWhirter^  and  D'Angelo"'’  Both  these  models  treat  re¬ 
combination  of  the  electron  into  a  bound  excis  .d  state  and  are  called  radia- 
tional  recombination.  The  treatmen*  -eference  6  for  hvdrogen-like  ions  is 
the  most  precise  available  and  \ .  ill  be  used  in  this  report  both  as  v  asis  for 
flow-field  calculations  to  predict  the  electron  density  distribution  and  for 
comparison  with  experimental  rate  constants.  In  a  recent  report  Byron, 
Stabler,  Bortz®  showed  that  the  results  of  reference  6  could  be  obtained  readily 
with  some  simplifying  assumptions.  This  work  however  did  not  alter  the  phy¬ 
sical  model  or  numerical  values  of  reference  6. 

Two  recent  recombination  experiments  in  helium  plasma  have  been  carried 
cur  by  Hinnov  and  Hirschberg^  and  Robben,  Kunkel  and  Talbot^.  Both  groups 
of  authors  compared  their  data  with  reference  6.  The  agreement  was  excell¬ 
ent  in  both  papers. 

The  experiments  discussed  in  this  paper  treat  electron-ion  recombination  in 
a  thermal  cesium  olasma  over  a  temperature  range  from  4300  to  1200  degrees 
helium  at  electron  densities  from  1.  7  (1014)  to  10  ^  electrons  per  cubic  cen¬ 
timeter.  The  plasma  investigated  was  a  turbulent  nonequiiibr iu  a  wake  gen¬ 
erated  by  the  passage  of  a  hypersonic  free-flight  sphere  through  an  argon- 
cesium  vapor  mixture.  Electron  density  was  measured  by  two  microwave 
focused  probes  at  frequencies  of  117  and  35  kilomegacycles .  In  general,  the 
measured  recombination  rates  in  these  experiments  are  in  reasonable  agree¬ 
ment  with  the  theoretical  values  of  reference  6. 

In  this  report,  the  experimental  technique  and  equipment  will  be  described. 

This  includes  a  discussion  of  the  microwave  plasma  diagnostic  apparatus. 

The  f'.ow-field  features  as  well  as  computing  model  employed  for  determining 
temperature  and  evaluated  nonequilibrium  electron  density  distribution  are 
presented.  Finally,  the  electron  density  measurements,  recombination  rates, 
and  compar  son  with  theoretical  values  from  reference  6  are  given. 
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TABLE  J 


PREVIOUS  MEASUREMENTS 
OF 

RECOMBINATION  RATES  a  IN  A  CESIUM  PLASMA 


II.  EXPERIMENTAL  TECHNIQUE 


The  cesium  plasma  investigated  in  this  study  was  generated  thermally  behind 
the  bow  shock  wave  of  a  5.  67 -mm  diameter  hyper  velocity  steel  sphere  tra¬ 
velling  in  an  argon-cesium  vapor.  The  bore-sized  projectile  was  accelerated 
to  13.  000  ft/sec  by  compressed  hydrogen  in  a  light  gas  launcher  The  experi¬ 
mental  facility  is  sketched  schematically  in  figure  1.  A  detailed  description 
of  the  launcher,  range  apparatus,  and  experimental  piocedure  is  contained  in 
appendix  A 

The  projectile  travels  through  a  sit  all  tank  which  absorbs  the  blast  irom  the 
expanding  hydrogen  driver  gas.  Next,  it  parses  through  two  stations  with 
light  screens  for  measurement  of  the  projectile  velocity  and  generation  of 
electronic  pulses  for  triggering  the  sweeps  of  the  oscilloscope  for  the  down¬ 
stream  microwave  apparatus  The  blast  tank  is  filled  with  commercial  grade 
argon  at  the  same  pressure  as  used  in  the  test  section  Next,  the  projectile 
pierces  a  thin  stainless  steel  liaphram  (.0005  inches  thick)  and  commences  to 
traverse  the  10  feet  long  l-i/2-inch  diameter  pyrex  test  section.  This  heated 
test  section  is  previously  filled  with  reagent  grade  argon  and  chemically  pure 
cesium  vapor.  The  details  of  the  construction  and  operation  of  the  oven  unit 
as  well  as  the  gas  loading  cycle  are  contained  in  appendix  A.  Each  firing  re¬ 
quires  several  days  The  cesium  is  allowed  to  coat-out  the  interior  walls  of 
the  reaction  vessel.  This  method  was  established  by  Haught  (reference  il) 
using  a  heated  shock  tube  to  measure  onset  of  ionization.  Cesium  is  admitted 
to  the  test  section  is  liquid  form.  The  pressure  is  determined  from  the  vapor 
pres  sure -temperature  data  in  reference  (12). 

The  projectile  then  passes  in  front  of  an  unfiltered  radiometer  (photomultiplier 
type)  which  monitors  the  variation  with  time  of  the  level  of  visible  radiation 
from  the  projectile  and  wake  at  that  station.  This  is  a  qualitative  device  to 
establish  the  absence  ol  high-velocity  particles  which  might  arise  from  up¬ 
stream  events  (launch  tube  fragments,  diaphram  pieces,  etc  ).  Further  down¬ 
stream,  the  spherical  model  goes  through  the  focused  spot  (15-mm  diameter) 
of  the  )  17  kilomegacycle  (F  -band)  microwave  station.  After  another  foot  of 
travel,  the  projectile  arrives  at  the  focused  spot  of  the  35  kilomegacycle 
(K-band)  microwave  plasma  diagnostic  station.  Electron  density  is  measured 
at  both  stations. 

Two  additional  measurements  are  made  at  the  end  of  the  test  section.  These 
are  the  rotating  mirror  streak  photographs  and  the  spectra.  The  first  of  these 
provides  a  photographic  measurement  of  the  projectile  velocity,  establishes 
the  absence  (or  presence)  of  debris  in  the  wake,  and  indicates  the  extent  of 
wake  radiation  (which  results  showed  to  be  limited  to  a  region  within  three 
diameters  of  the  sphere).  The  spectra  are  used  to  check  the  purity  of  the 
gas  samples  used  in  the  program.  They  showed  only  lines  from  cesium  and 
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Figure  1  SKETCH  OF  THE  CHEMICAL  KINETICS  BALLISTIC  RANGE  FACILITY 


argon  with  the  sole  exception  of  the  sodium  resonance  line.  Rotating  mirror 
streak  and  spectra  measurements  are  discussed  in  appendix  B.  Finally,  the 
projectile  pierces  the  stainless  steel  exit  diaphram  and  is  ultimately  disinte¬ 
grated  in  the  steel  plates  of  the  arresting  device.  The  theoretical  model 
employed  for  analysis  of  the  microwave  measurements  and  the  microwave 
apparatus,  are  described  in  detail  in  appendix  C,  These  measurements  were 
based  on  the  low  refractive  index  n,  (o<n<  1.  0)  of  the  highly  ionized  cylindrical 
plasma  Using  a  ray-tracing  technique,  it  is  snown  in  appendix  C  that  the 
power  of  the  signal  received  at  the  detecting  horn  is  related  to  the  refractive 
index  (or  equivilent  electron  density)  in  the  plasma. 


III.  EXPERIMENTAL  MEASUREMENTS 


The  experimental  program  was  carried  out  in  a  gaseous  mixture  of  argon  and 
cesium  This  was  necessitated  by  the  fact  that  the  highest  usable  oven  temp¬ 
erature  was  limited  to  260° C  bv  the  volatility  of  several  teflon  seals.  This 
factor  limited  the  maximum  cesium  vspor  pressure  to  0.  7  iorr.  Calculations 
based  on  shock -tube  measurements  by  Haught^  on  the  rate  of  ionization  in 
argon-cesium  mixtures  indicated  that  for  5,  67-mm  diameter  spheres  at  0.  7 
torr  the  incubation  time  for  ionization  in  the  shock  layer  would  exceed  the 
residence  time  in  the  shock  layer.  Thus,  it  was  necessary  to  dilute  the  cesium 
with  an  inert  gas  in  ordei  to  increase  the  frequency  of  collisions  in  the  shock 
layer. 

The  danger  in  the  procedure  was  that  (a)  the  argon  might  ionize  and  (b)  with 
the  higher  argon  density  three -body  recombination  between  a  cesium  ion,  an 
electron  and  an  Argon  atom  might  dominate  the  electron  removal  process. 

Preliminary  experiments  were  carried  out  in  pure  argon  (16  electron  volts 
ionization  potential)  (reference  13)  at  50  torr.  These  demonstrated  negligible 
wake  ionization,  Several  higher  velocity  experiments  were  run  subsequently. 
Again  electron  density  in  pure  argon  at  50  torr  was  observed  to  be  negligible. 
However,  experiments  at  500  torr  produced  a  significant  ionization  contri¬ 
bution  from  the  argon.  It  was  concluded  that  the  experiments  could  be  run  at 
velocities  up  to  15,  000  ft/ sec  in  argon  at  pressures  up  to  50  torr  without 
contributions  from  the  argon  ionization. 

Discharge  tube  measurements  of  electron-ion  recombination  in  argon  and  neon 
plasmas  have  been  made  by  Kentyf14)  at  3000° K,  while  Biondi  and  Brownf  *5), 
and  Redfield  and  HoltOM  (17)  have  reported  measurements  at  300° K.  Their 
quoted  recombination  rates  fitted,  to  within  a  factor  of  three,  a  three -body 
ion-electron-atom  recombination  rate  of  (3)  (10"22)  cc2/(particle)2/sec 

T  "°"K 

Under  the  conditions  of  the  experiments  carried  out  in  the  investigation  re¬ 
ported  here,  this  three -body  mechanism  would  result  in  a  second-order  rate 
of  2(  1 0 ~ ® )  cc/elect.  /sec.  This  rate  is  at  least  two  order  of  magnitude  faster 
than  the  radiative  collision  rates  of  reference  (6)  for  cesium  plasmas.  Thus, 
the  experimental  measurements  would  readily  demonstrate  the  effect  of  the 
Argon  atoms  on  i ecombination  rate.  This  effect  will  be  discussed  further  in 
section  V. 

The  1 1 7-kilomegacycle  microwave  probe  detects  plasma  in  the  range  10l  3  < 
ne-  <  2  (10*4)  cc'1.  The  35  kilomegacycle  probe  operates  in  the  range 
1 0 1 2  <  ng  -  <  1 0 1  3 .  These  probes  are  separated  by  approximately  one  foot 
along  the  test  section.  The  data  from  each  probe  is  plotted  on  the  same  figure 
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for  each  run  This  provides  a  contiguous  record  of  the  electron  density-time 
history  in  the  range  cf  electron  densities  from  2(10^)  through  10^  Cc '  -1  All 

the  microwave  data  were  reduced  according  to  the  calibration  plots  derived  in 
appendix  C  Approximately  twenty  runs  were  made,  hut  many  had  to  be  rejected 
because  the  projectile  trajectory  deviated  by  more  than  a  projectile  radius  from 
the  focused  spot  of  the  microwave  probe,  or  spurious  particles  were  present 
in  the  wake  On  four  runs  the  projectile  impacted  and  damaged  the  test  section. 

Shown  in  figure  2  is  a  composite  of  the  experimental  data  for  the  0.  1  torr 
cesium  shots  The  vertical  bars  on  the  experimertal  graphs  indicate  the  un¬ 
certainty  in  measurements  These  uncertainties  are  discussed  in  appendix 
C  This  same  figure  also  shows  calculated  electron  density  distributions 
These  were  determined  usin^  the  nonequilibrium  computing  program  described 
in  appendix  D.  Two  calculated  curves  are  shown  which  correspond  to  upper 
and  lower  limits  in  the  initial  enthalpy  assumption.  This  assumption  is  dis¬ 
cussed  in  appendix  D. 

It  is  interesting  to  note  that  the  10-percent  variation  in  initial  temperature  in¬ 
herent  in  the  upper  and  lower  calculated  limits  did  not  cause  significant  devi  ’ - 
tion  in  the  electron  density  curves  The  calculated  and  experimental  curves 
are  very  similar  in  shape  but  displaced  by  about  a  factor  of  two  in  absolute 
magnitude  of  the  electron  density. 

Similar  results  were  obtained  for  the  0.  6  torr  cesium  data  shown  in  figure  3. 
Here,  the  calculated  electron  density  is  lower  by  approximately  a  factor  of 
three  from  the  data. 

Since  the  computed  electron  density  distribution  is  based  on  the  theoretical 
recombination  rates  of  reference  (6),  this  agreement  between  experiment  and 
calculation  is  indicative  of  agreement  between  experimental  recombination 
rates  and  the  theoretical  values  in  reference  (6).  In  the  next  section,  the 
electron  density  data  will  be  analyzed  to  determine  recombination  rates  and 
these  will  be  compared  with  values  from  reference  (6). 
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Figure  2  VOLUME  WAKE  ELECTRON  DENSITY  AS  A  FUNCTION  OF  TIME 
(P„.  ~  0.10  TORR) 


IV.  DETERMINATION  OF  REACTION  RATES 


A.  FLOW  FIELD 

In  the  chemical  kinetic  experiments  being  performed  in  this  program,  a  small 
amount  of  cesium  mixed  with  argon  is  heated  to  temperatures  greater  than 
15,  000°K  in  the  shock  layer  of  a  hypersonic  spherical  projectile.  In  equilib¬ 
rium,  the  cesium  would  be  completely  ionized,  while  only  10 '4  parts  of  the 
argon  would  be  ionized.  Thus,  initially  the  cesium  ions  will  dominate  the 
ionized  species.  It  may  be  expected  that  the  cesium  ions  will  dominate  the 
nonequiiibrium  regions  alt  .  Equilibrium  ionization  does  not  occur  instan¬ 
taneously  in  the  sho  :k  layer. 

Several  reactions  may  occur  as  the  gas  traverses  the  shock  layer.  These  are: 


A  +  A  -»  A+  +  e~  -t-  A  R(  1 ) 

A  4  e~  ■*  A+  +  e~  +  e~  R(2  ) 

Cs  +  A  *  Cs+  +  e“  +  A  R<3) 

Cs  -t  Cs  Cs+  +  e-  +  Cs  R(4) 

Cs  +■  e~  -*  Cs*  +  e~  +  R(5) 


Here  A  and  A+  represent  an  argon  molecule  and  ion;  Cs  and  Cs+  indicate  a 
cesium  molecule  ana  ion;  ie  the  symbol  for  an  electron. 

Reactions  R(l),  F  (3),  and  R  (4)  refer  to  atom-atom  collisions  which  result 
in  the  ionization  of  one  of  the  atoms.  The  cross  sections  for  these  reactions 
are  small  compared  with  those  for  reactions  R.(2)  and  R  (5)  resulting  from 
electron-atom  collisions.  However,  since  the  atom  concentrations  are  very 
large,  ionization  by  atom-atom  collisions  dominate  initially  (except  for 
reaction  R(4)  which,  due  to  the  small  cesium  concentration,  may  be  neglected). 
When  the  electron  concentration  is  within  a  few  orders  of  magnitude  of  equilib¬ 
rium,  electron-atom  processes  R(  2)  and  R  (5)  become  dominant. 

The  preliminary  experiments,  reference  (13),  were  carried  out  in  velocities 
lower  than  9500  ft/ sec.  Low  levels  of  wake  ionization  (109  cc“^)  were  ob¬ 
served.  This  was  attributed  to  failure  of  the  cesium  to  achieve  equilibrium 
ionization  in  the  shock  layer,  This  was  overcome,  in  the  set  of  experiments 
reported  here,  by  increasing  the  sphere  velocity  to  13,000  ft/s ec.  This  in¬ 
creased  the  shock  layer  temperature,  and  the  wake  electron  density  was  ob¬ 
served  to  be  in  excess  of  2  ( 10 ^)  cc ~ 3 . 
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In  the  wake  region,  mixing  of  the  high  temperature  core  with  the  cool  outer 
fluid  causes  the  temperature  to  decrease.  Recombination  of  electrons,  pri¬ 
me  rily  by  the  inverse  of  reaction  (Rfj),  is  the  mechanism  by  which  the  plasma 
adjusts  to  the  lower  temperature. 

Inasmuch  as  the  recombination  reactions  take  place  in  a  turbulent  wake,  the 
reduction  and  interpretation  of  the  experimental  data  require  a  theoretical 
model  of  such  a  nonequilibrium  wake. 

The  wake  produced  oy  a  ballistic  range  model,  moving  at  hypersonic  speed 
through  a  fluid  medium,  is  considered  to  be  composed  of  inviscid  and  viscous 
portions.  The  outer  wake  is  a  region  which  spreads  from  the  outer  edge  of 
the  inner  wake  across  the  wake  shock  to  the  bow  shock,  and  the  inner  wake 
is  a  strongly  decelerated  core  portion  of  the  flow  behind  the  body.  The  viscous 
inner  core  initially  becomes  very  hot  because  of  the  viscous  heating  in  the 
boundary  layer  over  the  body  surface  and  the  recompression  at  the  wake  neck. 

The  flow  passing  through  the  strong  portion  of  the  bow  shock  and  coming  into 
the  outer  inviscid  region  is  also  heated  to  high  temperatures.  Under  this 
condition,  ionization  initially  occurs  in  the  vicinity  of  the  model,  and  in  the 
downstream  (cooler)  region  of  the  wake,  recombination  will  begin.  It  is  the 
purpose  of  the  present  study  to  utilize  the  results  of  a  turbulent  wake -flow 
field  analysis  to  interpret  ballistics  range  experiments  concerning  the  recombi¬ 
nation  rate  in  a  cesium -argon  plasma. 

To  minimize  the  efforts  ind  the  cost  of  the  operation,  and  also  to  eliminate  the 
effect  of  ablation  as  well  as  the  effect  of  the  angle  of  attack,  spheres  havt  been 
used  as  the  testing  models  in  the  ballistic  range  experiments  fo*  the  present 
study. 

Unlike  the  wake  behind  a  slender  body,  where  the  boundary  conditions  along 
the  edge  of  the  inner  wake  are  uniform,  the  analysis  of  the  wake  behind  a  blunt 
body  requires  a  careful  treatment  of  the  ewailowing  of  the  hot  inviscid  outer 
wake  by  the  viscous  inner  wake.  Moreover,  the  pressure  gradient  terms  appear¬ 
ing  in  the  momentum  and  energy  equations,  which  are  usually  neglected  for  a. 
slender  body  case,  must  now  be  retained.  The  detailed  analysis  is  described 
in  appendix  D. 

B.  EVALUATION  OF  ENVIRONMENTAL  DATA 

The  agreement  between  calculated  and  measured  magnitudes  of  the  electron- 
density  decay  with  time  is  reasonable,  as  seen  from  figures  2  and  3.  How¬ 
ever,  the  recombination  rate  is  proportional  to  the  rate  of  change  of  the  elec¬ 
tron  density  with  time.  In  this  regard,  measurement  and  calculation  differ 
markedly.  Thus,  it  has  been  necessary  to  develop  a  procedure  for  determin¬ 
ation  of  the  rate  constants  directly  from  the  data.  It  is  desirable  to  fit  the  data 
with  a  second  order  rate  constant,  since  the  experimental  values  will  be  com¬ 
pared  with  theoretical  values,  reference  (6),  which  are  tabulated  in  second- 
order  form. 
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The  rate  of  change  of  electron  density  in  a  plasma  with  only  electrcn-ion 
recombination  processes  occuring  is 


2 


where  a  is  the  second-order  recombination  rate  in  cc/electron/  sec 


(0 


The  effect  of  diffusion  uas  evaluated  separately  by  calculating  several  frozen 
fiow  cases  using  the  computing  program  described  in  appendix  D.  These 
calculations  showed  that  diffusion  produced  a  negligible  electron  density  de¬ 
cay  in  the  time  interval  of  these  measurements  (t  *>  10"^  sec.  ).  Thus,  it  is 
satisfactory  to  account  for  the  electron  removal  by  only  recombination  pro¬ 
cesses. 


Integration  of  this  rate  equation  over  the  time  interval  from  t 


1 


to  t  yields 


a  (ne_  ,  T)  dt 


(2) 


For  a  sufficiently  small  time  interval,  the  recombination  may  be  considered 
constant.  Tien  equation  (2)  becomes 


(3) 


Values  of  a  were  determined  from  the  experimental  data  using  equation  (3) 
after  the  data  was  replotted  as  inverse  jetron  density  as  a  function  of  time. 
Such  a  plot  is  shown  in  figure  4.  Upper  and  lower  limits  for  the  rate  constants 
were  found  by  measuring  the  maximum  and  minimum  slopes  through  groups  of 


adjacent  data  points. 


An  initial  point 


c  required  to  evaluate  the 


first  measured  electron  density  data  point,  was  defined  by  an  equilibrium  wake 
flow  field  calculation.  This  was  taken  to  be  the  most  downstream  wake  loca¬ 
tion  for  which  the  cesium  could  remain  completely  ionized  (T,  =  5000°K).  Be¬ 
yond  tj.  the  calculated  temperature  and  associated  equilibrium  electron  density 
decay  rapidly  and  nonequilibrium  in  the  electron-ion  reactions  causes  the 
electron  density  to  persist  at  higher  than  equilibrium  levels.  It  is  only  in  this 

region  ‘hat  equation  (3)  is  valid.  The  calculated  initial  points  ar e(~j  =10 


cc/electron,  t.  =  200  microseconds  for  0.  1  torr  -  p 

i  cesium; 


1.  7  10 


-16 


cc/electron,  zOO  microseconds  for  0.  6  torr  =  ^cesium 
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The  rate  constants  reduced  vvith  this  method  are  shown  in  figures  (5)  through 
(H)  alone  with  theoretical  values  from  reference  (6).  Upper-and  lower-limit 
c  irves  inchoate  the  uncertainly  in  rate  constant  caused  by  the  spread  in  meas¬ 
ured  electron  density. 

The  theoretical  values  of  a  shown  in  figures  5  through  8  were  obtained  from  the 
theoretical  •  alues  of  reference  (6)  (shown  in  figure  9),  oy  using  the  electron 
density  from  the  measurements  and  temperature  calculated  from  the  turbulent- 
v.  ake  computing  program  described  in  appendix  D,  The  calculated  temperature 
distribution  is  shown  in  figure  10.  Upper-and  lower-limit  temperatures  are 
indicated,  which  correspond  to  the  maximum  and  minimum  initial  enthalpy 
approximations  used  in  the  calculation.  This  approximation  is  discussed  in 
more  detail  in  appendix  D.  The  theoretical  values  of  a  presented  used  com¬ 
binations  of  the  largest  electron  density  (in  the  measured  uncertainty  band) 
with  the  lower -limit  temperature  and  the  lowest  electron  density  with  the 
upper-limit  temperature.  This  resulted  in  the  largest  uncertainty  band  in  the 
theoretical  valns  of  a  .  The  cross-hatched  region  on  figure  9  indicates  the 
temperature -electron  density  range  covered  in  this  investigation. 

In  general,  for  the  0.  1  torr  cesium  data  (figures  5-7)  the  experimental  values 
of  the  rate  constants  agree  reasonably  with  the  theoretical  recombination  rates. 
The  measured  rates  tend  to  be  faster  than  the  theoretical  by  about  a  factor  of 
two.  For  the  0.  fc  torr  cesium  data  (figure  8)  the  experimental  values  fall  with¬ 
in  the  tl  eoretical  band  in  the  low-temperature  range  (T<  1900*X). 

Comparison  of  experiment  with  theory  is  shown  in  figure  11  for  all  the  data. 
Here,  the  ratio  of  midrange  experimental  to  midrange  theoretical  rate  con¬ 
stant  ta'  er,  from  figures  5  through  8  is  indicated  as  a  function  of  temperature. 
There  does  not  appear  to  be  a  specific  trend  in  the  ratio  a  “Bates  with 

t‘"opera*  lire. 
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Figure  6  COMPARISON  OF  THEORETICALLY  AND  EXPERIMENTALLY  DETERMINED 
RECOMBINATION  COEFFICIENTS  FOR  ROUND  4-B 
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Figure  8  COMPAP  SON  OF  THEORETICALLY  AND  EXPERIMENTALLY  DETERMiNED 
RECOMBINATION  COEFFICIENTS  FOR  ROUND  4-D 
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EC0M6INATI0N  COEFFICIENT,  a  C  C  /,„ctron /seCOnd 


Figure  9  RECOMBINATION  COEFFICIENT  FROM  BATES,  KINGSTON,  McWHIRTER, 
(REFERENCE  6)  AS  A  FUNCTION  OF  TEMPERATURE 
AND  ELECTRON  DENSITY 
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Figure  10  CENTERLINE  WAKE  TEMPERATURE  AS  A  FUNCTION  OF  DISTANCE 
BEHIND  SPHERE;  VELOCITY  -  1300  FT/SF.C. 


4000 


VALUES  OF  RATE  CONSTANT  AS  A  FUNCTION  OF  TEMPERATURE 


V.  SUMMARY 


The  results  of  an  experimental  determination  of  electron-cesium  ion  recombin¬ 
ation  rate  have  been  presented.  The  cesium  was  ionized  thermally  by  the  bow 
shock  wave  of  a  spherical  proiectile  travelling  13,  000  feet  per  second.  In  the 
wake  turbulent  wave  behind  the  sphere,  the  plasma  decays,  since  the  temperature 
decreases  rapidly  and  recombination  becomes  important.  The  recombination 
rate  was  determined  from  the  measured  wake  electron  density  decrease  wi'h 
time.  F -band  and  Ka-band  focused  microwave  probes  were  used  in  this 
measurement.  The  electron  density  distribution  was  in  reasonable  agreement 
with  calculations  from  a  nonequilibrium  turbulent  wake  computing  program. 
Similarly,  the  recombination  rates  obtained  from  the  electron  density-time 
data  were  in  agreement  with  the  theoretical  values  of  Bat"**,  Kingston,  and 
McWhirter  based  on  collisional  radiative  recombination.  Values  of  the  re¬ 
combination  rate  constant  were  determined  ranging  from  1.  3  ( 10"  ^)cc/elect/ 
sec  at  a  temperature  and  electron  density  of  3400°K  and  1.  4  (10^)  elect 

- 1 0  1 2  CC 
respectively  to  (8.  0)  10  cc/elect/sec  at  1250eK  and  2  (10  )  elect/cc.  The 

rate  constant  was  approximately  twice  of  the  Bates  value  throughout.  The 

temperatures  used  were  calculated  values  from  the  wake  program. 

The  microwave  and  ballistic  range  equipment  and  techniques  used  in  these  ex¬ 
periments  are  discussed  in  deta  1  in  the  appendixes  to  this  report. 
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APPENDIX  A 


THE  CHEMICAL  KINETICS  BALLISTICS  RANGE 


The  experimental  investigation  has  been  carried  out  in  a  ballistic  range  designed 
tor  chemical  kinetic  studies.  As  such  it  contains  many  innovations  in  design  as 
compared  to  conventional  range  facilities,  and  its  operation  requires  applica¬ 
tion  of  new  techniques.  The  design  and  operation  of  this  facility  is  described 
in  the  appendix. 

1.  DESCRIPTION  OF  FACILITY 

The  facility  is  composed  of  a  light  gas  gun  for  Lunching  the  model,  a  blast  tank 
for  absorbing  the  gun  blast,  a  test  section,  and  »n  arresting  device  for  stopping 
the  projectile.  Auxiliary  equipment  includes  a  light  screen  for  measuring  the 
projectile  velocity  and  necessary  instrumentation  to  make  experimental  meas¬ 
urements.  The  instrumentation  and  its  use  is  described  in  some  detail  in 
appendixes  B  and  C. 

The  test  section  in  which  the  measurements  are  made  consists  of  a  pyrex  glass- 
pipe  lange  tube  in  which  the  atmosphere  and  pressure  are  controlled,  and  is 
surrounded  by  a  temperature -controlled  oven.  Additional  equipment  consists 
of  a  vacuum  system  for  exhausting  the  test  chamber  and  the  range,  a  gas  infec¬ 
tion  system  for  injecting  high-purity  gase  1  into  the  test  chamber,  and  an 
ampoule  crusher  for  introducing  high-purity  metal  vapors  into  the  test  chamber. 

A  photograph  of  the  entire  facility  is  shown  in  figure  12.  The  individual  com¬ 
ponents  are  described  in  detail  in  the  following  praagraphs. 

a.  Light-Gas  L.  icher ,  Blast  Tank,  and  Arresting  Device 

The  light-gas  launcher  (figure  13)  is  a  two -stage  ballistic  launching  device 
which  uses  gunpowder  to  drive  an  aluminum  piston  into  the  first-stage  pump 
tube  which  is  filled  with  a  light  gas  such  as  hydrogen  or  helium,  under  a 
pressure  of  500  psi.  The  piston  compresses  the  gas,  which  in  turn  ruptures 
a  shear  disc  at  the  entrance  to  a  second-stage  barrel  that  contains  a  bore- 
size  spherical  projectile.  The  gas  propels  the  projectile  through  the  barrel 
and  into  the  test  arge.  Velocities  up  to  15,000  fps  are  thus  attained, 
using  type  440  stainless  steel  spheres.  An  expendable  section  with  a  con¬ 
stricted  bore  is  employed  between  the  two  stages  to  trap  the  piston. 

The  circular  steel  four-petal  shear  disc  has  two  transverse  serrations 
which  begin  to  rupture  cleanly  at  their  intersection  and  then  fold  forward 
into  a  curved  conical  section  at  the  breech  of  the  second-stage  barrel. 

This  technique  was  developed  to  prevent  pieces  of  the  shear  disc  from  con¬ 
taminating  the  experiment.  The  projectile  is  held  by  friction  fit  at  the 
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Figure  12  PHOTOGRAPH  OF  THE  ENTIRE  FACILITY 


Figure  13  SKETCH  OF  LIGHT  GAS  LAUNCHER 


exit  of  the  conical  section,  until  it  is  dislodged  by  the  pressure  pulse  from 
the  ruptured  shear  disc. 

The  blast  (absorbing)  tank,  shown  in  figure  1,  is  a  six-foot  long  cylinder 
with  a  diameter  of  eighteen  inches. 

The  projectile  arresting  device  s  a  16-inch  length  of  8-inch  diameter  pipe 
coupled  to  the  end  of  the  test  chamber  by  glass  pipe.  It  contains  a  series 
of  increasing  thickness  steel  plates,  backed  by  1-inch  armor  plate.  A 
0.035-inch  aluminum  plate  is  placed  foremost  to  determine  the  point  of 
projectile  impact  to  be  used  as  a  check  on  the  projectile  trajectory  align¬ 
ment. 

b.  Light  Screens 

Two  light  screens,  spaced  three  feet  apart,  are  mounted  on  quad¬ 
ruple-port  aluminum  couplings.  They  are  of  the  dark-field  type, 
where  a  tungsten  light  source  is  projected  through  the  projectile  flight 
path,  to  a  low- reflectance  mat  on  the  opposite  side  of  the  range  test  section. 
A  photomultiplier  tube  is  directed  at  th’  mat  and  receives  a  lignt  pulse 
only  by  reflection  from  the  passing  projectile.  The  light  screen  signals 
are  used  to  start  and  stop  an  Events  Per  Unit  Time  counter  and  to  trigger 
the  operation  of  the  microwave  apparatus. 

c.  Environmental  Test  Section 


The  test  chamber  is  a  Pyrex  glass  tube  l-l/3-inch  diameter  by  108  inch 
long,  through  which  the  projectile  is  fired.  It  is  separated  from  the  r  st  of 
the  range  by  diaphragms  placed  at  each  end.  A  3/4-inch  by  16 -inch  verti¬ 
cal  glass  pipe  stem  is  located  near  one  end  of  the  Pyrex  pipe  to  permit 
evacuation  and  for  insertion  of  test  gases.  The  stem  is  coupled  to  the 
vacuum  system  by  a  glass  stopcock.  A  fixture  on  the  stem  provides  the 
means  for  crushing  a  glass  ampoule  in  vacuo,  thereby  permitting  the 
injection  of  suitable  materials  into  the  test  chamber.  Heating  of  the 
ampoule  crusher,  stopcock,  and  vertical  stem  by  electrical  heating  tapes 
is  necessary  to  prevent  cesium  condensation  in  these  parts  which  protrude 
from  the  oven. 

The  oven,  figure  (14),  is  a  130-inch  long  reinforced  case  with  a  23-inch 
square  cross-section,  and  is  made  of  1 /4-inch  sheet  aluminum.  It  encloses 
an  8-inch  square  cross  section  heating  chamber  made  of  1-inch  thick 
aluminum  plate,  (see  figure  15.)  The  apace  between  the  two  is  filled  with 
fiberglass  insulation.  The  heating  chamber  rests  on  6  inches  of  Marinite 
slab  insulation.  A  removable  6-inch  thick  insulated  cover  serves  to  close 
both  the  heating  chamber  and  case.  The  entire  oven  is  mounted  on  two 
screw  jacks  which  serve  to  raise  and  lower  it  (figures  16  and  17).  In  the 
lowered  position, work  may  proceed  unhampered  on  the  glass  test  chamber. 
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Figure  14  SKETCH  OF  OVEN 
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Figure  15  CROSS-SECTIONAL  SKETCH  OF  OVEN 


igure  16  PHOTOGRAPH  OF  OVEN  ''RAISED) 
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Figure  1  7  PHOTOGRAPH  OF  OVEN  (LOWERED' 


The  one -inch  thick  bottom  plate  of  the  heating  chamber  has  ZOChromalox 
heating  elements  spaced  evenly  along  its  length.  These  elements  are 
installed  in  holes  bored  transversely  through  the  plate.  Initial  heating  of 
the  elements  causes  them  to  expand  and  contact  tne  walls  of  their  noles, 
thereby  increasing  heat  conduction.  The  total  powei  consumption  of  the 
20  elements  is  7000  watts.  The  elements  are  paired  by  means  of  switches 
to  provide  10  heat  zones.  Each  pair  of  elements  has  three  heat  positions, 
HIGH,  LOW,  and  OFF,  corresponding  to  350  watts,  175  watts,  and  0  per 
ehment.  This  permits  the  maintenance  of  various  temperature  gradients 
along  the  heating  chamber.  Thermal  pumping  of  the  medium  within  the 
test  chamber  may  then  take  place.  Oven  temperature  may  be  controlled 
from  ambient  to  300°C.  Oven -tempe rature  protection  is  provided  by  a 
detector  imbedded  in  the  bottom  plate  and  another  inside  the  oven  proper. 
Heat-up  time  from  ambient  to  200°C  is  approximately  2  hours.  Tempera¬ 
ture  uniformity  (to  *  2  C‘)  is  established  in  an  additional  hour.  Along  the 
sides  of  the  oven  and  in  line  with  the  glass  test  chamber  are  four  pairs  of 
double  -  glas  sed  ports  and  one  single  port.  These  may  be  used  for  visual 
inspection  of  the  range  or  for  instrumenation.  The  inner  glass  ports  are 
of  Pyrsx,  The  outer  ports  are  selected  plate  glass.  Chromel -Alumel 
thermocouples,  a  reference  cold  junction,  and  a  Leeds  and  Northrup 
recorder  are  used  to  monitor  temperatures. 

The  vacuum  system  is  in  two  sections.  The  first  section  is  a  mechanical 
pump  for  evacuating  the  blast  tank,  pump  tube,  projectile  arresting  device, 
and  range  piping  exclusive  of  the  test  chamber.  The  second  section  used 
to  (evacuate)  the  test  chamber,  consists  of  a  mechanical  roughing  pump, 
three-stage  oil  diffusion  pump,  and  cold  traps.  The  test  chamber  is 
separated  from  the  rest  of  the  range  by  0.005-inch  stainless  steel  diaphrams 
just  inside  the  ends  oi  the  oven. 

Independent  operation  of  the  mechanical  pumps  permits  evacuation  of  the 
range  proper  and  test  chamber,  so  that  a  minimum  pressure  differential 
exists  across  the  frangile  diaphragms.  V iton  O  rings  are  used  as  inter¬ 
face  glass-to-glass  seals.  A  large  ground-glass  ball  joint  is  used  to  inter¬ 
connect  the  oil  diffusion  pump  and  the  liquid  Nitrogen  cold  trap.  All  vacuum 
components  are  readily  demountable  for  cleaning.  A  minimum  of  vacuum 
grease  is  used.  A  Wallace  and  Tiernan  absolute  pressure  gauge  is  used 
to  indicate  range  pressure.  A  Coons  diaphragm  vacuum  gauge  is  used  to 
indicate  rough  test  chamber  pressure.  A  Pirani  thermal  conductivity 
gauge  and  a  McLeod  gauge  are  used  to  indicate  test  chamber  pressure. 

A  thermocuple  vacuum  gauge  is  used  to  indicate  roughing  pump  pressure. 

The  ampoule  crusher  is  an  oyvaen-free  high-conductivity  copper  tube  v,ith 
a  screw  jack  at  one  end  and  a  range  coupling  adapter  at  the  other.  Glass 
heating  tape  is  used  to  elevate  its  temperature  A  glass  ampoule  of  the 
correct  material  is  inserted  in  the  crusher  when  assembling  the  range.  At 
the  appropriate  time,  the  screw  jack  is  activated,  thereby  crushing  the  glass 


-34- 


ampoule,  and  permitting  the  material  !t  contains  to  drop  down  into  the 
test  chamber. 

2.  EXPERIMENTAL  TECHNIQUE 

It  is  essential  that  the  projectile  je  fired  through  'he  renter  of  the  glass  test 
chamber  within  ±  3  mm.  To  do  this,  the  barrel  is  aligned  on  a  paper  target  by 
means  of  a  theodolite.  A  powder  breech  is  then  affixed  to  the  barrel  and  a  bore 
size  Zelux  projectile  is  fired  at  4000  fps  velocity.  Corrections  in  barrel 
orientation  are  made  until  it  is  shooting  centrally. 

The  barrel  is  then  cleaned  and  the  remainder  of  the  light  gas  gun  is  assembled 
to  the  barrel  without  changing  varrel  orientation.  This  condition  is  verified 
by  means  of  four  dial  indicators  which  are  in  contact  with  the  barrel.  Soft 
copper  washers  are  used  to  seal  all  interfaces  of  the  light  gas  gun. 

The  ends  of  the  Pyrex  chamber  are  fitted  with  stainless  steel  diaphragms  and 
Teflon  gaskets.  The  chamber  is  coupled  to  the  range  at  both  ends  with  spring- 
loaded  flanges,  which  allow  for  any  thermal  expansion  as  well  as  for  any  Teflon 
flow  at  elevated  temperatures.  The  range  piping  outside  the  oven  is  fitted  with 
Teflon  expansion  bellows  to  allow  lor  longitudinal  expansion  of  the  test  chamber, 
as  well  as,  minor  manufacturing  variations  in  tubing  length. 

Several  different  diaphragm  materials  were  tested  to  find  the  most  suitable  one. 
Some  materials  were  subject  to  premature  rupture  due  to  the  corrosive  Teflon 
sheet  overlays  and  Teflonized  aluminum  were  also  tr’ed.  These  materials  did 
not  corrode;  however,  the  presence  of  carbon  lines  in  the  slit  spectrograph 
indicated  that  contamination  by  vaporizing  Teflon  was  taking  place.  Ultimately, 
stainless  steel  diaphragms  of  0.0005-inch  thickness  were  employed. 

An  additional  contaminant  was  discovered  in  the  reaction  of  Silicone  stopcock 
grease  with  cesium.  A  stainless  steel  bellows  valve  with  Viton  seals  was  tried 
in  place  of  the  greased  stopcock.  Failure  of  the  Viton  when  maintained  a* 

200°C  for  prolonged  periods  forced  its  removal  from  the  system.  Ultimately, 
the  problem  was  alleviated  by  returning  to  the  glass  stopcock  and  carefully  wiping 
out  the  excess  grease  after  the  stop  cock  plug  had  been  inst.  fed  and  seated.  In 
the  low-temperature  glass  piping  external  to  the  oven,  it  was  found  that  Viton 
O  rings  formed  a  better  vacuum  seal  than  the  Teflon  gaskets  normally  used 
with  the  commerical  flanges. 

The  cesium  vapor  is  introduced  into  the  test  chamber  from  a  2.  5-gram  ampoule 
of  cesium.  The  ampoule  is  placed  in  the  crusher  and  sealed  with  a  piston, 
located  at  the  end  of  the  screw,  which  contains  two  Viton  O  rings.  The 
ampoule  must  be  free  of  all  labeling  substances  to  avoid  introducing  contaminant 
into  the  system. 
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A  rough  vacuurr  is  drawr.  on  the  test  chamber  and  range  simultaneously,  a. id 
the  system  is  checked  for  leaks.  A  minimum  pressure  differential  is  maintained 
across  the  diaphragms. 

When  the  test  chamber  pressure  is  100  microns  of  less,  the  three  -  stage  oil 
diffusion  pump  is  activated.  After  attaining  a  further  decrease  in  pressure 
to  b  microns  or  less,  the  liquid  nitrogen  traps  are  filled.  A  mass  spectrometer 
is  then  employed  to  detect  any  remaining  leaks  in  the  test  chamber.  A  satis¬ 
factory  condition  is  attained  when  the  leak  rate  is  i.ss  than  1  micron  per  minute, 
without  cryogenic  or  oil  diffusion  pumps  operating. 

The  lower  portion  of  the  oven  is  then  elevated  to  enclose  the  test  chamber,  and 
the  oven  lids  are  put  in  place.  Marinite  collars  are  used  to  support  the  test 
chamber  in  the  oven.  A  two-piece  cupola  is  used  to  enclose  the  vertical  range 
stem.  Glas s -asbestos  heating  tape  is  wound  about  the  ampoule  crusher 
assembly  and  range  stopcock.  The  oven  and  the  heating  tapes  are  turned  on. 
Temperatures  of  both  oven  and  heating  tape  are  indicated  by  Chromel-Alumel 
thermocouples . 

Reagent  grade  argon  is  injected  into  the  test  chamber  to  a  pressure  of  5  torr 
gauge  for  one  minute,  after  which  the  chamber  is  again  evacuated.  This  serves 
to  flush  the  chamber. 

After  a  final  check  of  the  leak  rate,  the  test  chamber  stopcock  is  closed  and 
the  ampoule  crus*-  r  activated.  The  cesium  drops  into  the  test  chamber, 
forming  a  small  pool  in  the  r  mge  pipe.  The  minimum  amount  of  cesium 
required  to  conduct  the  experiment  is  0.  01  grams.  The  ampoule  contains  2.  50 
grama,  thereby  assuring  a  sufficient  quantity  of  cesium  in  the  test  chamber. 

The  cesium  will  coat  the  inner  walls  of  the  test  chamber.  The  partial  pressure 
of  cesivm  is  not  measured  directly,  but  is  determined  from  the  vapor  pressure 
temperature  table  (reference  12),  the  temperature  being  measured  by  the  six 
Chromel-Alumel  thermocouples  located  along  the  test  section. 

Following  the  cer  um  injection,  and  at  periods  of  7  minutes,  20  minutes,  1  hour, 
and  at  1-hour  intervals  thereafter,  a  check  is  made  to  determine  the  impurity 
gas  generation  rate  caused  by  reaction  of  some  cesium  with  impurities  in  the 
system.  To  determint  the  gas  gene"ation  rate  at  each  of  these  times,  the  test 
section  is  evacuted.  The  hign  vacuum  system  stopcock  is  then  shu.  off  and  the 
test  chamber  stopcock  is  opened,  leaving  the  test  clumber  connected  to  the 
Piraru  gauge.  The  pressure  reading  is  noted,  and  a  series  of  readings  taken 
at  10-second  intervals  foi  1-1/2-minute  period.  This  series  of  readings, 
figure  18,  constitute  the  gas  generation  rate,  and  are  indicative  of  the  mpun- 
ties  remaining  in  the  system.  The  test  chamber  heating  cycle,  is  continued 
until  the  impurity  is  less  than  100  microns  at  the  two-minute  check  peiiod. 

The  vacuum  system  stopcock  is  closed,  sealing  off  the  test  chamber  and  gauges. 
Reagsnt  gmde  argon  is  injected  to  the  d<  sired  pressure.  Five  minutes  are 
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Figure  13  GAS  GENERATION  RATE 


allowed  for  the  argon  to  diffuse  throughout  the  chamber.  The  range  is  charged 
with  sufficient  air  to  reduce  the  pressure  differential  across  the  diaphragms  to 
less  than  1  torr  gauge. 

At  this  point  the  proper  g<<  s  mixture,  pressure,  and  temperature  have  been 
established  in  the  test  chamber.  The  projectile  may  now  be  fired  and  the  data 
recorded. 
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APPENDIX  D 


OPTICAL  MEASUREMENTS 

1.  ROTATING  MIRROR  STREAK  PHOTOGRAPH 

An  Avco  Mode'  MC  300-1  rotating  mirror  camera  is  used  to  record  the  possible 
presence  oi  debris,  such  as  diaphram  or  shear  disc  fragments  which  may  follow 
the  projectile  down  the  range.  The  appearance  of  any  such  fragments  would 
signify  the  presence  of  impurities  and  hence,  invalidate  the  data  collected  on 
the  test  firing.  Also,  the  camera  provides  a  check  on  the  projectile  velocity 
measurement  made  from  the  light  screens.  Figure  19  showns  a  debris  free  test 
firing,  while  figure  20  showns  a  test  firing  which  has  debris. 

In  operation  the  camera  receives  an  image  from  the  obje.  t  source  and  sweeps 
the  image  by  means  of  a  rotating  mirror  onto  fast  emulsion  film  strips.  The 
resulting  trace  presents  a  direct  position  versus  time  record  of  the  events 
being  studied.  This  camera  utilizes  a  magnetic  tachometer  to  accurately 
measure  the  rotational  speed  of  the  mirror,  which  in  turn  determines  the 
writing  speed. 

2.  OPTICAL  SPECTA 

Primary  interest  in  the  optical  spectral  was  the  application  cesium  line  broaden¬ 
ing  (Stark  effect'  to  measure  electron  density  optically.  Since  the  Stark  con¬ 
stants  for  cesium  lines  have  veen  evaluated  by  Criem #,  this  technique  offered  a 
useful  comparison  to  the  microwave  diagnostics.  This  theory  shows  that  even 
for  the  most  sensitive  lines,  line  widths  must  be  resolved  to  approximately  one 
Angstrom  for  determination  of  electron  densities  of  lC^.cm*^.  This  resolu¬ 
tion  was  obtainable  with  available  spectrographs. 

A  preliminary  test  was  conducted  to  see  if  the  intensity  of  the  cesium  spectra 
was  sufficient  to  obtain  spectra  with  an  f6  time  -  resolved  Tropel  spectrograph. 
The  results  of  these  measurements  showed  that  the  intensity  was  not.  sufficient 
to  obtain  an  exposure. 

Spectral  data  were  obtained  with  a  transmission  grating  spectrograph,  however, 
this  spectrograph  is  non-time  resolving;  i.  e.  .  it  integrates  all  the  radiation 
from  the  projectile  and  its  wake.  The  increased  exposure  levels  obtained 
through  the  use  of  a  wide  slit  (500  microns)  allowed  spectral  data  to  be  recorded 
at  the  expense  of  spectral  resolution.  Table  II  shows  the  relative  intensities  of 
spectral  lines  for  four  test  firings  in  the  facility.  A  considerable  enhancement 
of  the  cesium  spectrum  is  shown  when  the  cesium  pressure  is  increased  from 
0.0  to  0.7  torr.  Figure  21  shows  some  of  the  spectra  taken  by  the  spectrograph. 
While  these  spectra  were  interesting,  they  were  inadequate  as  a  backup  electron 

density  determination. 

•  " 

Grtem,  H.  R  ,  P!  >m.i  Spectroscopy .  McGraw-Hill,  New  York 
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Figur«  19  ROTATING  MIRROR  STREAK  PHOTOGRAPH  -  CLEAN  SHOT 
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Fi<jur«  20  ROTATING  MIRROR  STREAK  PHOTOGRAPH  -  SHOT  WITH  DEBRIS 
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I  ABLE  II 


WAVELENGTHS  AND  RELATIVE  INTENSITIES  OE  SPECTRAL  LINES 


X 

W’a  velength 

1  (Angstroms) 

Element 

Cs  =  Cesiun 

A  =  Argon 
Na  =  Sodium 

Relative 

ntensity 

for 

50  torr  A 

0.  1  torr  Cs 

Relative 

Intensity 

for 

500  torr  A 

0.  1  torr  Cs 

Relative 

Intensity 

for 

50  torr  A 

0  torr  Cs 

i 

Relative 

In'e  asil\ 
for 

50  torr  A 

0.  7  f err  Cs 

|  3876.  39 

Cs 

F 

’  ' 

3888.  65 

Cs 

F 

3947.  5 

A 

F 

)948.  98 

A 

F 

4044.  42 

A 

1 

F 

F 

4158.  59 

A 

F 

F 

F 

F 

4164.  18 

A 

F 

F 

F 

F 

4181.  88 

A 

I 

I 

F 

F 

j  4191.03 

A 

I 

1 

F 

F 

4198.  32 

A 

I 

1 

F 

E 

j  4200.67 

A 

I 

I 

F 

F 

4251.  18 

A 

I 

I 

1 

IF 

4259. 36 

A 

I 

T 

I 

IF 

4266.  29 

A 

T 

I 

I 

IF 

4272.  17 

A 

I 

I 

F 

IF 

4333.  56 

A 

F 

F 

F 

F 

4335. 34 

A 

F 

F 

F 

F 

4555. 35 

Cs 

IB 

4593.  18 

CS 

I 

5162 

A 

I 

I 

I 

I 

5187 

A 

I 

I 

I 

I 

5257 

Cs 

F 

5340 

Cs 

I 

5407 

Cs 

I 

5414 

Cs 

I 

5466 

Cs 

IB 

5503 

Cs 

F 

IB 

556  3 

Cs 

F 

F 

5566. 7 

Cs 

F 

F 

56  35 

Cs 

F 

B 

4 

Cs 

F 

B 

5746 

Cs 

F 

5844. 7 

Cs 

F 

F 

B 

588  7.  6  2 

A 

F 

F 

F 

|  5888. 59 

A 

F 

F 

F 

5890 

Na 

F 

B 

5896 

Na 

F 

B 

6010 

Cs 

I 

F 

B 

6032.  12 

A 

F 

F 

JT 

6034. 09 

Cs 

F 

F 

B 

6043.  23 

A 

F 

F 

F 

oToZTTZ 

A 

F 

F 

F 

6059.  37 

A 

F 

F 

F 

6212.  9 

Cs 

I 

I 

B 

6  355 

Cs 

F 

F 

_ 1 

Code:  B  -  bright,  I  ~  intermediate,  F  =  faint,  IB  =  intermediate  bright, 
IF  =  intermediate  faint. 
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Figure  21  TIME  INTEGRATED  SPECTRA 
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APPENDIX  C 


MICROWAVE  MEASUREMENTS 


The  use  of  microwave  techniques  for  plasma  diagnostic  studies  is  common  in 
hypervelocity  ballistic  ranges.,  (see,  references  1-3).  Although  there  are 
various  methods  for  measuring  the  plasma  properties,  attention  here  will  be 
confined  to  the  microwave  free-space  focused  probe.  The  focused  probe  can 
be  used  to  determine  the  magnitude  and  spatial  distribution  of  wake  ionization 
characterized  by  the  electron  density  and  collision  frequency  of  electrons  with 
neutral  atoms  and  ions.  These  later  quantities  can  be  determined  for  mam' 
simple  plasma  configurations  by  measuring  the  transmission  coefficients  of  the 
plasma.  Transmission  coefficients  are  the  absorption  coefficient  a,  and  the 
phase  coefficient  /9  .  These  constants,  listed  below,  have  been  derived  from  the 
equation  of  motion  of  an  electron  under  the  influence  of  a  sinusoidal  driving 
force  (reference  4). 


where 

-  plasma  frequency  =  -  eVm<o 

ne  =  electron  density 


v  =  collision  frequency 

oj  =  signal  frequency 


c  =  charge  on  an  electron 


m  =  mass  of  an  electron 


c  =  velocity  of  light  in  free  space 


=  permittivity  of  free  space 
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tor  the  gas  densities  employed  in  these  experiments,  the  calculated  collision 
frequency  '’■eterem.es  5  and  b)  is  negligibly  small  compared  'o  signal  frequencies 
of  15  to  117  kmc  (see  figure  22).  Therefore,  equations  1  and  2  reduce  to 


For  any  plasma  diagnostic  apparatus,  the  usable  signal  frequency  range  is 
primarily  governed  by  the  plasma  frequency  cl,  of  the  plasma  being  examined, 
since  for  practical  measurements,  the  signal  frequency  cannot  differ  too  greatl’ 
from  the  pla  sma  frequency.  This  is  apparent  from  the  interaction  equations 
(3)  and  (4)  above,  viz,  when  the  plasma  frequency  is  equal  to  the  signal  fre- 

Cdp 

quenc^  or,  -  =  1,  then  a  =  =  0;  and  when  the  plasma  frequency  is  very 

0) 

CJ 

small  a>  «  (jj,  a  =  0,  ft  =  — 

P  c 

To  measure  electron  densities  of  10^  to  10^  elect/cc,  as  encountc  red  in  the 
cesium  recombination  experiments,  the  microwave  signal  frequency  must  be 
in  the  35-  to  1  1 7 -kilorr.egacy cle  range. 

Since  the  collision  frequency  does  not  appear  in  equation  4,  the  electron  density 
may  be  solved  for  from  only  the  phase  coefficient  a.  General  Motors,  *  RCA  ^ 
and  RARDE^  have  applied  this  method  to  obtain  electron  density  through  the 
use  of  a  microwave  focused  probe  interferometer.  The  interferometer  meas¬ 
ures  the  shift  in  phase  of  the  transmitted  signal  in  the  plasma.  From  this  phase 
shift  the  phase  coefficient  is  obtained  and  the  electron  density  is  calculated.  To 
be  able  to  use  this  type  of  measuring  system  the  following  conditions  must  be 
met: 

1.  The  microwavo  energy  must  be  concentrated  so  that  only  a  negligible 
part  of  the  beam  may  bypass  the  plasma. 

2.  In  the  immediate  area  of  the  focal  plane,  ’.he  constant  phase  contours 
are  planes  which  are  approximately  parallel  to  the  focal  plane. 

3.  The  electromagnetic  wave  front  must  be  normal  to  a  plane,  parallel 
plasma  slab.  (Plasma  diameter  must  be  much  greater  than  the  projectile 
diameter). 

These  conditions  would  not  all  be  satisfied  in  the  present  test  facility,  because 
with  the  5.  6b -mm  diameter  spheres,  the  cylindrical  wake  plasma  does  not  meet 
criteria  (2)  or  (3).  Criterion  (1)  is  occasionally  violated  since  the  normal 
projectile  trajectory  deviations  are  of  the  order  the  projectile  radius. 
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Several  types  ot  microwave  apparatus  could  be  employed  lor  the  necessary 
measurements.  F  or  example,  transverse  parallel  focused  microwave  probes 
discussed  in  reference  7  would  overcome  the  difficulty  of  the  projectile  tra¬ 
jectory  deviation.  Unfortunately,  microwave  components  were  not  available 
tor  this  type  apparatus.  Consequently,  a  different  microwave  attenuation 
measuring  technique  was  developed  for  installation  in  the  chemical  kinetics 
range.  In  this  method, measurements  of  electron  density  from  ?(10^)  to  10^ 
elect/cc  are  made  using  two  focused  free-space  probes  operating  at  frequencies 
of  35  kmc  and  117  kmc.  The  probes  are  aligned  normal  to  the  projectile  tra¬ 
jectory.  Each  probe  consists  of  a  klystron  transmitter,  transmitting  horn  antenna, 
tv  o  Teflon  dielectric  focussing  lenses,  receiving  horn  and  crystal  detector  (see 
figure  23).  One  lens  focuses  the  transmitted  energy  to  the  axial  center  of  the 
ballistic  range  and  the  other  lens  collects  and  refocusses  the  beam  to  the 
receiving  horn  antenna.  The  focus  spot  size  of  the  beam  for  both  probes  is 
16  mm  iti  diameter.  This  is  large  enough  to  encompass  any  deviations  form 
the  projectiles  normal  flight  path,  since  the  wake  diameter  in  the  region 

X 

( — - >  150) of  measurement  is  greater  than  28  mm  (i.  e.  ,  nearly  twice  the  spot 

d 

diamete  r) . 

The  probes  measure  the  attenuation  in  the  transmitted  signal  due  to  the  com¬ 
bined  effects  of  the  absorption  coefficient  a,  the  phase  coefficient  ft.  (see 
Equations  3  and  4)  This  measured  attenuation  is  a  function  of: 

1.  Attenuation  in  the  plasma  due  to  the  absorption  coefficient. 

2.  Refraction  or  defocussing  of  the  focused  beam  due  to  the  phase 

coefficient. 

3.  Reflection  from  the  plasma  boundaries  due  to  the  phase  coeff  cient. 

By  applying  the  calculated  theoretical  collision  frequencies  (figure  22)  and 
their  corresponding  electron  densities  to  equation  3,  the  attenuation  caused  by 
can  be  calculated.  The  attenuation  in  the  plasma  as  a  result  of  is  shewn  in 
figure  24  as  a  function  of  electron  density  for  signal  frequencies  of  both  35 
anc.  177  kmc.  The  index  of  refraction,  n  ,  (  n=  c  ft  /  <v)  attenuates  the  transmitted 
signal  by  causing  refraction  at  the  plasma  boundaries. 

Index  of  refraction  versus  electron  density  for  33-  and  117-kmc  signal  fre¬ 
quencies  is  shown  in  figure  25.  Assuming  a  homogeneous  plasma,  Snell's  law 
of  refraction  is  applicable.  Also,  if  the  plasma  is  assumed  to  h  ve  a  known 
cylindrical  shape,  geometric  ray  tracing  can  be  applied  (reference  8).  As  the 
electron  density  increases,  the  index  of  refraction  decreases  towards  0,  and  the 
number  of  rays  w-hich  are  deflected  away  from  the  receiving  horn  increases. 

This  effect  is  illustrated  schematically  in  figure  26.  By  applying  the  rav-trac- 
ing  technique  to  plasma  cylinders  of  varying  indices  of  refraction,  the  attenua- 
t’-n  due  to  refraction  as  a  furction  of  electron  density  is  obtained  (see  figure  27). 
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Figure  23  SCHEMATIC  DIAGRAM  OF  FOCUSED  FREE  SPACE  MICROWAVE 
PROBE  INSTALLATION 


ATTENUATION  IN  DB  FOR  12.8-MM  LENGTH  OF  PLASMA  AS  A  FUNCTION 
OF  ELECTRON  DENSITY 


65-9387 
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Figure  25  INDEX  OF  REFRACTION  AS  A  FUNCTION  OF  FLECTPON  DENSITY 
AT  35  AND  117  KMC 


Figure  26  ILLUSTRATION  OF  THE  EFFECT  OF  MICROWAVE  REFRACTION  BY 

PLASMA  CYLINDER 
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Figure  27  ATTENUATION  BY  REFRACTION  OF  RAYS  BY  CYLINDRICAL  PLASMA 


Reflection  at  the  plasma  boundar.es  causes  attenutation  in  the  transmitted 
signal  also.  The  reflection  coefficient  as  a  function  of  electron  density  for  35 
and  117  kmc  is  shown  in  figure  28,  where  the  reflection  coefficient  is  given  by 

Reflection  coefficient  = 

The  total  attenuation  is  obtained  from  the  triple  product  of  the  attenuations 
caused  by  the  absorption  coefficient,  refractive  index,  and  the  reflection 
coefficient.  The  curve  of  this  triple  product  is  shown  in  figures  28  and  29  for 
the  two  signal  frequencies. 

To  obtain  the  greatest  accuracy  from  this  measuring  system,  electron  density 
data  points  were  usually  only  taken  in  the  regions  where  the  index  of  refraction 
was  nearly  zero  and  unity.  This  procedure  yields  four  electron-density  data 
points  per  test  firing;  two  from  the  35-kmc  probe  at  2  x  10^  and  1.  5  x  10^ 
el  -  cc  “  ^  and  two  from  the  117  -  kmc  probe  at  1  x  10  * 3  and  1.7  10^el-cc_^. 

The  inaccuracies  in  this  system  due  to  the  assumptions  that  the  plasma  is 
homogeneous,  isotropic,  and  of  known  size  have  been  taken  into  account  in  the 
electron  density  profile  curves  by  presenting  the  data  points  and  their  inaccur¬ 
acies  as  error  bands.  In  addition  to  uncertainties  in  the  data  from  noise  level 
in  the  oscillograms,  certain  errors  are  inherent  in  the  microwave  techniques. 

The  error  bands  were  estimated  in  the  following  way: 

1.  The  35-kmc  and  the  117-kmc  probe  data  points  were  compared  to  each 
ciher  in  the  region  of  overlapping  sensitivity. 

2.  Three-kmc  microwave  cavity  plasma  measurements  from  earlier  pure 
argon  results  shots  were  compared  to  the  35-kmc  and  the  117-kmc  probe 
results. 

3.  Errors  were  determined  in  measurements  of  dielectric  rods  of  known 
refractive  index. 

The  limits  on  the  error  bands  estimated  by  this  procedure  are  as  follows: 

117 -kmc  probe - n  «■  0  - 1.3  to  1.7^  el cc*  3 

n  «.  1 - 0.  8  to  2.  013  el  cc*3 

35-krrc  pro'  - n  0 - 1.  1  to  1.  5^3  el  cc-3 

n  «,  1 - 1.0  to  4.0^  el  cc-3 
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Figure  28  REFLECTION  COEFFICIENT  AND  ATTENUATION  FROM  REFLECTION  AS  A 

FUNCTION  OF  ELECTRON  DENSITY 
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The  following  is  a  listing  of  the  components  used  in  the  microwave  circuits 
(figure  23) 


35-kmc  Equipment  177-kmc  Equipment 

Micro  power  supply  mod.  FXR  mod.  Z315B 
dfd-cs-em-av 


K q ui pment  Description 
Klystron  power  supply 

Klystron 

Isolation 

Horn  antenna 

Lenses 

Crystal  detector 
Oscilloscope 


VarionAssoc.  V0239A 
400  mw  -  35  kmc 

Cascade  research 
Mod,  KA-13-1 

DeMornay  Bonardi 
DBD  520  10  db  horn 

Avco-made  Teflon 
lenses 

Microwave  Assoc. 
Mod.  539c 

Tektronix  Mod.  502 


VarionAssoc.  V0715D 
140  mw  -  117  kmc 

TRG  mod.  110-41 

DeMornay  Bonardi 
DBW  520  15  db  horn 

Avco-made  Teflon  lenses 

Microwave  Assoc. 
MA4147 

Tektronix  Mod.  502 


To  assure  stability  of  the  tramsmitter,  the  klystron  was  isolated  from  the  rest 
of  the  circuit  by  a  ferrite  isolator. 

The  focusing  lenses  are  made  of  Teflon  to  withstand  the  cesium  oven  tempera¬ 
tures.  Surfaces  of  the  lenses  are  hyperbolic  to  avoid  spherical  aberrations. 
The  lenses  are  of  nonsymmetrical  design  with  the  surfaces  facing  the  horns 
naving  a  focal  length  of  3.  0  inches  and  the  surface  facing  the  plasma  having 
a  focal  length  of  3.  5  inches.  The  field  distribution,  or  focus  spot  size,  was 
measured  by  passing  two  identical  small  spheres  through  the  focal  plane  for 
various  spacings  between  the  spheres.  Where  the  two  spheres  are  separated 
by  more  than  the  beamwidth,  the  field  distribution  in  the  focal  plane  is  traced 
out  by  each  sphere.  As  the  spheres  are  placed  closer  together,  the  two  pat¬ 
terns  ultimately  overlap.  The  focused  spot  size  is  defined  as  the  separation 
distance  between  the  two  identical  spheres  in  the  focal  plane  which  produces 
field  patterns  that  overlap  that  overlap  at  the  3-db  points.  The  focal  spot  size 
for  the  lenses  used  measured  16  mm  in  diameter. 


A  precision  calibrated  attenuator  is  used  to  calibrate  the  output  voltage  of  the 
crystal  detector  receiver  in  terms  of  relative  input  power.  The  output  voltage 
is  displayed  vertically  on  an  oscilloscope,  while  the  scope's  horizontal  trace 
is  swept  in  time,  resulting  in  a  display  of  attenuation  versus  time.  Since  the 
velocity  of  the  projectile  is  known  from  the  light  screen  measurements,  the 
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time  can  be  converted  into  distance  and  a  plot  of  attenuation  versus  distance 
behind  the  projectile  can  be  obtained.  The  sweep  of  the  oscilloscope  is  triggered 
externally  from  the  light  screens  via  a  time  delay  trigger  generator.  The  delay 
of  the  trigger  generator  is  adjusted  to  allow  the  projectile  to  travel  from  the 
light  screen  station  to  a  position  a  few  inches  ahead  of  the  microwave  probe 
stations  before  the  scope  is  triggered.  This  results  in  the  best  possible  time 
resolution  by  allowing  faster  scope  sweep  rates  to  be  used.  The  oscilloscope 
trace  is  recorded  with  a  poloroid  camera. 

Electron-density  profiles  are  obtained  by  applying  the  calibration  curves, 

(figures  29  and  30),  to  the  oscilloscope  attenuation  record.  This  procedure  is 
illustrated  in  figure  31.  As  previously  mentioned,  the  electron  density  data 
points  are  taken  in  the  region  where  the  index  of  refraction  is  approximately 
0  and  1,  for  it  is  in  this  region  that  the  greatest  accuracy  is  obtained.  Inter¬ 
mediate  points  are  taken  in  those  few  cases  where  a  small  particle  followed 
too  closely  behind  the  sphere  to  permit  the  attenuation  to  return  to  the  zero-db 
level.  The  riples  in  the  output  signal  as  the  plasma  decays  is  the  result  of 
phase  shift  in  the  microwave  signal  passing  through  the  cylindrical- shaped 
plasma.  The  portion  of  the  signal  that  passes  nearer  the  edges  of  the  plasma 
through  the  center  of  the  plasma  cylinder.  Therefore,  the  microwave  signal 
emerging  from  the  plasma  consists  of  several  phase  fronts.  How  these  phase 
fronts  combine  vectorially  after  leving  the  plasma  will  determine  the  attenuation 
deviation  from  the  calibration  curves.  Sine'  these  attenuation  perturbations 
are  small,  they  have  only  a  minor  effect  on  the  accuracy  of  the  measuring 
system  and,  hence,  are  neglected. 
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Figure  29  Kg-BAND  PROBE  CALIBRATION 
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Figure  31  CORRESPONDENCE  BETWEEN  MICROWAVE  MEASUREMENTS  AND  THE 

REDUCED  DATA 
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APPENDIX  D 

THEORETICAL  TURBULENT  WAKE  ANALYSIS  WITH  ELECTRON-ION 

RECOMBINATION 


i 

i 

i 


I 
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Inasmuch  as  the  recombination  reactions  take  place  in  a  turbulent  wake  where 
the  coupling  of  chemical  an  .  iluid  dynamic  effects  is  significant,  the  interpreta¬ 
tion  of  the  experimental  daca  requires  a  theoretical  solution  of  such  a  non- 
quilibrium  wake.  In  this  appendix,  a  physical  model  is  developed  and  the 
method  of  solution  is  described. 


The  wake  produced  by  a  ballistics  range  model,  moving  at  hypersonic  speed 
through  a  fluid  medium,  is  considered  to  be  composed  of  inviscid  and  viscous 
portions,  as  illustrated  in  figure  32.  The  outer  wake  is  a  region  which  spreads 
from  the  outer  edge  of  the  inner  wake  across  tl  vake  shock  to  the  bow  shock, 
and  the  inner  wake  is  a  strongly  decelerated  cc  p:  portion  of  the  flow  behind 
the  body.  The  viscous  inner  core  initially  be  cc  Ties  very  hot  because  of  the 
viscous  heating  in  the  boundary  layer  over  the  body  surface  and  the  recom¬ 
pression  at  the  wake  neck.  The  flow  passing  through  the  strong  portion  of 
the  bow  shock  and  coming  into  the  outer  inviscid  region  is  also  heated  to  high 
temperatures.  Under  this  condition,  ionization  initially  occurs  in  the  vicinity 
of  the  model,  and  in  the  downstream  (cooler)  region  of  the  wake  recombination 
will  begin.  It  is  the  purpose  of  the  present  study  to  utilize  the  results  of  a 
turbulent  wake-fiow  field  analysis  to  interpret  ballistics  range  experiments 
concerning  the  recombination  rate  in  a  cesium-argon  plasma. 

To  minimize  the  efforts  and  the  cost  of  the  operation,  and  also  to  eliminate 
the  effect  of  ablation  as  well  as  the  effect  of  the  angle  of  attack,  spheres  have 
been  used  as  the  testing  models  in  the  ballistic  range  experiments  for  the 
present  study.  It  must  be  pointed  out  here  that  the  analysis  of  the  wake  be¬ 
hind  a  blund  body  is  not  exactly  the  same  as  the  analysis  of  the  slender  body 
case.  The  effects  of  the  swallowing  and  the  pressure  gradient,  which  are 
usually  neglected  for  a  slender  body  case,  may  not  be  neglected  for  a  blunt 
body.  However,  under  certain  conditions  of  the  present  experiments,  these 
two  effects  are  considered  unimportant.  This  conclusion  is  made  based  upon 
the  following  observations: 


(1)  The  theoretical  study  of  Lees  and  Hromas^,  indicates  that  the  process 
of  swallowing  enthalpy  or  momentum  defect  originally  contained  in  the 
outer  wake  is  almost  completed  in  about  300  diamete  s.  From  the 
Schlieren  photographs  of  Reference  2  the  swallowing  under  the  present 
experimental  conditions  is  found  to  be  completed  in  the  vicinity  of  150 
diameters. 


(2)  It  is  known,  from  characteristic  calculations,  that  the  effect  of 
pressure  gradient  seldom  extends  beyond  the  region  of  100  diameters. 
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Specifically,  the  distance  to  the  point  where  the  pressur*  in  the  wake 
returns  to  ambient  is  approximately  25  diameters  for  M  =  8.  5  and 
50  diameters  for  M  ^  =  22. 

(3)  The  fully  developed  turbulent  far -wake  is  usually  preceded  by  a  short 
region  of  laminar  wake.  The  transition  location  varies  according  to  the 
local  Reynolds  and  Mach  numbers.  Over  the  range  of  the  present  ex¬ 
periments,  transition  is  found  to  occur  prior  to  the  completion  of  the 
swallowing. 

(4)  In  the  prssent  analysis,  the  theoretical  computation  starts  from  the 
location  where  the  swallowing  is  completed.  The  initial  conditions  for 
the  turbulent  wake  development  are  matched  to  those  obtained  from  the 
experiments.  The  detailed  discussion  on  matching  conditions  is  given 
later.  To  start  the  theoretical  computation  from  the  completion  of 
swallowing  is  necessary  not  only  for  the  simplification  of  the  theoretical 
treatment,  but  also  for  matching  the  experimental  conditions. 

The  electron  density  prior  to  the  completion  of  swallowing  is  much  higher 
than  the  range  of  the  instrumentation  used  in  the  present  experiments. 

Thus  ,  it  is  corcluded  that  the  effects  of  swallowing,  pressure  gradient,  and 
the  laminar  wake  development  are  all  contained  in  the  portion  of  the  wake  prior 
to  the  initial  point  of  the  experiment,  and  therefore  they  are  excluded  in  the 
turbulent  wake  analysis  which  is  described  in  the  following  paragraphs. 

The  basic  equations  for  the  flow  of  a  chemically  reacting  multicomponent  gas 
mixture  are: 

<J(pj)  i  <?(pvr)  /iv 

- - +  —  -  =  0  '  ' 

&x  r  <Jr 
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The  notation  ueed  in  these  equations  is  contained  in  the  list  of  symbols.  The 
governing  equations  lor  turbulent  flow  are  obtained  when  concentration,  velocity, 
and  enthalpy  are  decomposed  into  their  mean  and  fluctuating  parts.  After 
time  averging,  they  reduce  to  a  set  of  equations 


where  is  the  turbulent  eddy  diffusivity  and  Prc  is  the  turbulent  Prandtl  num¬ 
ber.  For  simplicity,  the  usual  barred  notation  denoting  averaged  quantity 
(e.  g.  ,  u)  has  been  dropped.  In  deriving  the  above  enuations.  it  is  assumed 
that  the  turbulent  diffusion  is  predominant  over  the  laminar  diffusion  and  the 
value  of  turbulent  Lewis  number  is  unity.  According  to  the  measurements 
made  by  Fage  and  Falkner  in  the  wake  of  a  cylinder,  and  Reichardt  in  a  free 
jet,  both  of  which  are  described  in  reference  3,  the  magnitude  of  a  turbulent 
Prandil  number,  Prc  ,  was  found  to  vary  from  0.  5  to  0.  8.  A  value  of  0.  7  is 
used  throughout  the  present  computations. 

It  is  well  known  that  the  treatment  of  the  turbulent  problem  is  necessarily 
semiempirical.  The  knowledge  of  the  variation  of  some  physical  quantities, 
such  as  the  eddy  viscosity  and  the  correlation  functions,  is  primarily  obtained 
from  experimental  data.  The  magnitude  of.  eddy  viscosity  in  the  prese.L 
analysis  is  estimated  from  a  combination  of  mixing -length  theory  and  the  ex¬ 
perimental  results  of  Townsend^  as 


tp  -  c.b.  ( Au) 

where  C  is  a  proportionality  constant,  defined  as  the  reciprocal  of  the  universal 
eddy  Reynolds  number,  b  is  the  half  width  of  the  wake,  and  u  is  the  mean  vel  - 
ocity  defect  across  the  wake. 

Instead  of  trying  to  solve  the  governing  equations  exactly,  a  modified  integral 
method  is  applied.  In  terms  of  the  displacement  thickness  (  5*)  the  mo¬ 
mentum  thickness  (  )  ),  and  the  enthalpy  transfer  thickness  (  S  h)  defined  as 
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the  momentum  integral  and  the  energy  relations  become 


de 

(11) 

_  _  Q 

dx 

dSh 

(12) 

_ _  q 

dx 

which  can  be  solved  provided  that  suitable  forms  of  the  velocity  and  enthalpy 
profiles  are  assumed.  In  the  present  approach,  polynomials  of  fourth  degree 
are  chosen,  i.  e.  , 

u  —  u 

2  3  4 

t  -  -  =  a  +  a,  7)  +  a2  tj  +  J  +  a*  77’ 

ue  ~  uo 

and 

h  -  he 

g  55  T - Z~  ”  eo  +  el  ?  +  e2  I2  +  e3  ^  +  e4  ^ 

ho  *  "e 

Eight  of  the  10  free  constants,  a  ,  a, . a.te  ,  e. . e, 

o  1  4  o  1  4 

can  be  determined  from  the  following  eight  boundary  conditions. 

v 

At  if  -  0:  f  -  0,  g  -  1. 

t)  -  1:  f  ■  1,  g  r  0, 

The  velocity  and  enthalpy  profiles,  therefore,  can  be  written  as 
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(13) 


t  »  Arp  -  Irp  +  A  (17^  -  2  ip  +  ip) 

&  -  1  -  Arp  +  Irp  ¥  a{rp  -  2rp  +  rp)  (14) 


where  A  and  a  are  usually  called  the  shape  factors.  These  parameters.  A(x) 
and  a  (x)  ,  are  determined  by  satisfying  the  momentum  and  energy  equations 

along  the  wake  axis  in  the  forms 


u 


o 


duo 

17" 


and 


u 


o 


dho 

dx 


2  (, 


d2u 

77 


0 


(h„  -  h0)  -1- 
b2 


(15) 


(16) 


with  the  subscript  o  indicating  the  value  at  the  wake  axis.  The  integration  of 
species  conservation  equations  across  the  wake  can  be  treated  in  a  similar  way. 
They  can  be  reduced  to 


^  due 

u,  dx  dx 


— ~  rj dr;  (17) 
P 


where  the  concentration  exchange  thickness  for  species  i  is  defined  by  the  re¬ 
lation 


(ce  -  (Ce).  1  ydy 


The  concentration  profiles  for  species  i  are  derived,  similar  to  the  velocity 
and  enthalpy  profiles,  as 


1-4  p  +  3  rp  +  {rp  -  2  rp  +  p) 


(18) 


with  the  shape  parameters  ^  determined  by  the  species  continuity  equation 
along  the  wake  axis 
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09) 


*'v  d  C. 


Pf  dy'" 


<V. 


4,i-  f  -1 

-  —  [Voh  -  (CcM  - 

Pr‘  L  J  b2 


b*  +  P, 


<V. 


The  mass  rate  of  production  of  species  i  per  unit  volume  by  chemical  reactions, 
W.,  which  appears  everywhere  in  the  species  continuity  equation,  is  deter¬ 
mined  from  the  Law  of  Mass  Action 


V; 


M; 


dA; 
d  t 


-  (vf  -  *-)kf  n  (Aj)  '  +  (v:  -  vpkb  n  (Aj) 1 


(20) 


i  =*  l 


i  =  l 


where  A.  i3  the  molar  concentration  (mole/cc)  and  v',  v"  ,  are  the  stoichio¬ 
metric  coefficients  for  the  reaction  products  and  reactants,  respectively. 

The  complete  wake  solution  can  be  obtained  by  comgining  equations  (11)  through 
(19)  together  with  an  expression  for  species  production  (or  deletion)  from 
equation  (20). 

The  chemical  reactions  being  studied  in  the  ballistics  range  at  Avcc  RAD  in¬ 
volve  the  ion-electron  recombination  of  the  pure  alkali  metals.  In  the  first 
series  of  such  experiments,  5.67-mm  diameter  sphere  is  fired  into  argon  gas 
mixed  with  cesium  vapor.  Since  the  argon  cannot  react  with  the  cesium,  the 
only  reactions  which  will  be  observed  after  expansion  into  the  wake  will  be 

k, 


Cs+  +  e  +  (e  ) 


Cs  +  (e  ) 


II  Cs+  +  e“  +  (A) 

‘III 


Cs  +  (A) 


^In 

III  CsT  +  e—  - ♦  Cs  +  hi' 


Cs+  +  e" 

Applying  equation  (20)  for  the  reactions  of  (21)  yield? 


(21) 


and 


dA 

~d 


P-k.[v]  [v] 


(22) 


dA 

~d 


f~ =  -  )ki  K-] +  k”  M  +  km  j  [Ac-] 


K  +  k 


ii  M  j  [v] 


(231 
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where 


K*  “  kl  Ac-  +  kIII 


(24) 


The  combined  rate  constant  K*,  which  may  be  called  the  collisional  radiative 
recombination  coefficient,  has  been  determined  from  a  statistical  theory  in 
reference  5.  Its  value,  tabulated  as  a  function  of  temperature  and  electron  den¬ 
sity,  is  reproduced  figure  9.  As  to  the  three-body  recombination  rate  K  only 
limited  information  is  available  at  the  present  time.  However,  the  importance 
of  argon  as  a  third  body  can  be  tested  experimentally  by  varying  the  concentra¬ 
tion  in  the  ballistics  range  at  fixed  velocity  and  cesium  concentration.  From 
the  expression  (23),  it  can  be  seen  that  the  electron  density  decay  rate  is  con¬ 
trolled  by  argon  density,  if  K*«K^(A^).  Under  the  experimental  conditions 
of  the  present  study,  the  argon  density  is  kept  relatively  low  and  the  electron - 
density  decay  rate  is  found  to  be  little  affected  by  the  argon  density.  Therefore, 
theterm  K^(A^)  is  neglected  in  comparison  with  K*,  and  the  theoretical  elec¬ 
tron  distributions  in  the  wake  are  made  with 


K*  [Ae_  ]2 


(25) 


where  the  va’ues  of  K*  are  obtained  form  figure  (9). 

Computational  Procedure 


Following  the  method  of  theoretical  analysis  presented  above,  a  computer  pro¬ 
gram  was  written  co  facilitate  thenumerical  integration  by  IBM  7094  computer. 
The  initial  conditions  for  the  computation  are  specified  by  the  experiments  as 
follows: 

(1)  The  computation  starts  from  the  point  where  the  swallowing  io  com¬ 
pleted.  This  is  determined  from  the  Schiieren  photographs  of  reference  2. 

(2)  The  initial  electron  density  is  matched  to  the  experimental  value. 

(3)  The  wake  growth  rate  is  also  matched  to  the  experimental  value  meas¬ 
ured  from  the  Schiieren  photographs  of  reference  2. 

The  only  initial  condition  which  involves  some  uncertainty  is  the  initial  enthalpy. 
Since  measured  temperature  data  are  not  available,  the  initial  value  is  obtained 
from  the  approximated  calculation  of  reference  6.  To  include  the  effect  of  un¬ 
certainty  in  the  initial  temperature  estimation,  two  rune  representing  upper 
and  lower  limits  are  computed  for  ea  :h  case.  The  possible  range  of  the  solu¬ 
tion  is  thereby  determined. 
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Figures  2  and  3  in  the  text  respectively  indicate  the  center  line  electron  density 
in  the  wake  with  the  cesium  pressure  of  0.  1  and  0.  6  torr.  The  comparison  of 
these  two  figures  reveals  that  the  variation  of  initial  cesium  affects  mostly  the 
upstream  portion  of  the  wake,  while  the  downstream  electron  concentrations 
are  essentially  the  same.  This  behavior  is  explained  as  follows:  The  higher 
initial  cesium  concentration  (for  0.  6  torr),  resulis  in  a  highe~  initial  electron 
density.  From  figure  9,  it  is  seen  that  the  recombination  rate  K*  is  also 

higher  and  hence,  from  equation  (25),  -----  is  increased.  Thus  the  original 

d  t 

difference  in  electron  density  in  the  vicinity  of  the  wake  neck  »s  largely  elimin¬ 
ated  downstream  by  the  faster  initial  recombination.  Also  shown  in  figure  2 
and  3  the  effects  of  initial  temperature.  Based  on  previous  experience,  a 
±5  percent  uncertainty  may  be  expected  in  the  calculated  values  of  initial 
temperatu :e.  It  can  be  seen  that  this  difference  induces  only  small  differences 
in  electron  density.  This  insensitivity  of  wake  electron-aensitv  to  the  initial 
temperature  occurs  because  a  lower  initial  tempera*ure  results  in  a  higher 
recombiration  rate.  This  causes  the  electron  density  to  decrease  faster  than 
it  would  for  a  higher  initial  temperature  (see  figure  8)  Hence,  only  small 
differences  in  recombination  rates  and  electron  densities  are  ultimately  pro¬ 
duced. 

The  variation  of  the  temperature  along  the  wake  axis  is  presented  in  figure  9. 
The  temperature  distribution  is  practically  unaffected  by  the  cesium  pressure 
because  the  concentration  of  cesium  is  much  smaller  than  the  argon  in  the 
present  study.  Hence,  the  results  for  both  pressures  plot  as  single  curves. 
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